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Abstract Based on the most complete list of the results of an individual comparison of 
the proper motions for stars of various programs common to the Hipparcos catalog, each 
of which is an independent realization of the inertial reference frame with regard to stellar 
proper motions, we redetermined the vector lo of residual rotation of the ICRS system 
relative to the extragalactic reference frame. The equatorial components of this vector were 
found to be the following: u> x = +0.04 ± 0.15 mas yr' 1 , u y = +0.18 ± 0.12 mas yr -1 , and 
uj z = —0.35 ± 0.09 mas yr' 1 . 
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INTRODUCTION 

A kinematic analysis performed previously (Bobylev 2004) revealed appreciable residual 
rotation of the stars of the Hipparcos catalog (ESA 1997) with respect to the inertial refer¬ 
ence frame around the Galactic y axis with an angular velocity of —0.36 ± 0.09 mas yr' 1 
(milliarcseconds per year). A shortcoming of the method used is the absence of a rigorous 
criterion for separating the actual rotation of stars close to the Sun from the sought sys¬ 
tematic rotation. The method of analyzing the proper-motion differences for stars of various 
programs common to the Hipparcos catalog, each of which is an independent realization of 
the inertial reference frame with regard to stellar proper motions, is free from this short¬ 
coming. We call this method astrometric. Based on this method, Kovalevsky et al. (1997) 
found that all three components of the vector u> of residual rotation of the Hipparcos catalog 
with respect to the extragalactic reference frame have no significant deviations from zero, 
with the error of this vector along the three axes being ±0.25 mas yr' 1 . The difficulties in 
applying this method stem from the fact that there are few VLBI observable radio stars, 
while photographic catalogs are not free from the magnitude dependence of the stellar proper 
motions (the magnitude equation). The goal of this study is to redetermine the vector lo 
using the most complete list of the independent results from an individual comparison of the 
stellar proper motions. 

CHARACTERISTIC OF INDIVIDUAL SOLUTIONS 

Based on the NPM photographic program of the Lick observatory (Klemola et al. 1987), 
two catalogs of absolute proper motions for northern-sky stars have been published: the 
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NPM1 catalog of stars in 899 areas outside the Milky Way zone (Klemola et al. 1994) and 
the NPM2 catalog of stars in 347 areas of the Milky Way zone (Hanson et al. 2003). The 
areas of the NPM1 and NPM2 catalogs do not overlap. The results of a comparison of 
the proper motions for NPM1 and Hipparcos stars were described by Platais et al. (1998b) 
and designated as NPM(Yale) by Kovalevsky et al. (1997). An independent analysis of the 
proper-motion differences between the NPM1 and Hipparcos catalogs was also performed in 
Heidelberg; Kovalevsky et al. (1997) designated the results of this analysis as NPM (Heidel¬ 
berg). All of the above authors agree that the component u> z is difficult to determine using 
the NPM1 catalog, because there is a magnitude equation in the NPM1 proper motions that 
causes the NPMl-Hipparcos differences to be shifted by ~ 6 mas yr -1 for the brightest 
(~ 8 m ) and faintest (~ 12 m ) stars. Kovalevsky et al. (1997) did not include the component 
u z determined from NPM1 stars in their final solution. 

Zhu (2003) compared the proper motions of NPM2 and Hipparcos stars. This author 
found no appreciable magnitude dependence of the NPM2- Hipparcos differences. 

Based on the SPM program of photographic observations of southern-skv stars (Platais 
et al. 1995), Platais et al. (1998a) published the SPM2 catalog of absolute proper motions of 
stars in 156 southern sky areas. The results of a comparison of the absolute proper motions 
of SPM stars in 63 areas with the Hipparcos catalog were described by Platais et al. (1998a). 
We designate these SPM stars as SPM1. Zhu (2001) compared the proper motions of SPM2 
and Hipparcos stars. 

Based on the combined photographic catalog GPM (Rybka and Yatsenko 1997a), Rybka 
and Yatsenko (1997b) published a list of proper motions for bright stars common to the 
Hipparcos catalog designated as GPM1. Kislyuk et al. (1997) compared the proper motions 
of GPAI1 and Hipparcos stars. These authors concluded that there is a magnitude equation in 
the GPM1 catalog, and the parameters u x , u y , and u ; 2 determined only from faint comparison 
stars must be used. 

Bobylev et al. (2004) determined the parameters u x , uj y: and ou z by comparing the Pulkovo 
photographic catalog PUL2 and Hipparcos and found no appreciable magnitude equation in 
the proper motions of PUL2 stars. Each of the GPM and PUL2 catalogs is an independent 
realization of the photographic plan by Deutch (1954); the centers of the areas coincide and 
correspond to the list of Deutch (1955). 

The results of the well-known independent programs were presented by Kovalevsky et 
al. (1997). In addition, the results of the Potsdam program were presented by Hirte et al. 
(1996), and the Bonn program by Geffert et al. (1997) and Tucholke et al. (1997); the 
Earth Rotation Parameters (EOP) were analyzed by Vondrak et al. (1997); radio stars were 
analyzed by Lestrade et al. (1995), and the Hubble Space Telescope (HST) observations 
were analyzed by Kovalevsky et al. (1997). 

The main difference between our approach and the approach of Kovalevsky et al. (1997) 
is that we included the following data: 

(1) The results of a comparison of the proper motions for PUL2 and Hipparcos stars 
(Bobylev et al. 2004); 

(2) The results of a comparison of the proper motions for SPAI2 and Hipparcos stars 
(Zhu 2001); 

(3) The results of a comparison of the proper motions for NPM2 and Hipparcos stars 
(Zhu 2003); 

(4) The results of a comparison of the VLA absolute proper motions for radio stars and 
their Hipparcos proper motions. 
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Table 1 gives the components of the vector u> that we calculated using the absolute 
proper motions of 15 radio stars from Boboltz et al. (2003). These authors described in 
detail the project and the method of referencing the observations of radio stars to quasars 
and compared these stars with the Hipparcos catalog, but did not determine the components 
of the vector u). We used the equations in the form that was proposed and used by Lindegren 
and Kovalevsky (1995): 


A/j, a cos 8 — lj x cos a sin S + (1) 

+oj y sin a sin 5 — uj z cos 6, 

A fj,s = —c j x sin a + u y cos a , (2) 

where the Hipparcos catalog differences are on the left-hand sides of the equations. Here, 
<T 0 is the error per unit weight in the solution of Eqs. (1) and (2). The VLA (Boboltz et al. 
2003) and VLBI (Lestrade et al. 1995, 1999) absolute proper motions are independent. 

Our approach also differs from that of Kovalevsky et al. (1997) in that we used the 
results of a comparison of the NPM1 and Hipparcos stellar proper motions performed by the 
Heidelberg team, which were obtained from the data for 2616 stars, to determine uj z . The 
random errors in ui x , u y , and u z for the Heidelberg solution are given only for the sample of 
1135 stars (e Wa . = 0.25 mas yr _1 and e Uy = e Wz = 0.2 mas yr -1 ). The corresponding errors for 

the sample of 2616 stars are smaller by a factor of ~ ^2616/1135. Since there are problems 
with the NPM1 catalog, we assumed the random errors to be 0.25 mas yr -1 for u x and 0.20 
mas yr -1 for to y and uj z in order not to overestimate this solution. Our choice of this solution 
was dictated by the fact that it was obtained in the magnitude range 10.5 m — 11.5 m . As can 
be seen in Fig. lc from Platais et al. (1998b), the Hipparcos-NPMl proper-motion differences 
have a horizontal pattern near zero precisely in this magnitude range. In our opinion, the 
NPM1 proper motions are least affected by the magnitude equation in this magnitude range. 
In addition, we used only one of the available (not independent) results of a comparison of 
the SPM1 and Hipparcos catalogs, the YVA solution (Kovalevsky et al. 1997; Platais et al. 
1998b). 

DETERMINING THE VECTOR w 

We assigned a weight inversely proportional to the square of the error e in the corre¬ 
sponding quantities u x , u y . and uj z to each comparison catalog, which was calculated using 
the formula 

Pi = e kiev 2 /ei 2 , i = 1,..., 12, (3) 

where i is the number of individual sources, and e k i e v, the random error of the GPM1 (Kiev) 
program. Table 2 contains the data used here. The second column of the table gives the 
weights calculated using formula (3). Not all of the authors use equations in the form (1) 
and (2). In such cases, we reduced the signs of the quoted quantities (Zhu 2001, 2003) to 
the same form. 

The lower part of Table 2 gives the weighted means of u> x , tv y , and u> z . Mean 1 was 
calculated using all of the available data in Table 2. The SPM1 and SPM2 catalogs are not 
independent, but the methods of obtaining the solutions differ. The SPM1 solution (Platais 
et al. 1998b) was obtained only from the differences in /j a cos5 (e.g., only from Eq. (1)). 
The SPM2 solution (Zhu 2001) was obtained by simultaneously solving Eqs. (1) and (2); 
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the author pointed out that there is a color equation in the SPM2 catalog, which was not 
eliminated, and it is more pronounced in the (is differences. We calculated Mean 2 without 
using the SPM1 solution in order to analyze only independent sources. In this case, the 
SPM proper motions also have the largest weight when calculating the component u z , as 
can be seen from the table. A comparison of Means 1 and 2 reveals no significant differences 
between these two solutions, with lv z being an appreciable component. The errors of the 
vector <jj along the three axes e w = \Je Ux 2 ± e^ 2 + e^ 2 are ±0.21 and ±0.25 mas yr -1 
for Mean 1 and 2 solutions, respectively. Our solutions for the components u> x and u y do 
not differ significantly from the final solutions of Kovalevsky et al. (1997) based on both 
Linclegren’s and Kovalevsky’s methods. At the same time, there is a significant difference in 
the determination of u z , which we found to differ significantly from zero. The figure shows 
the projections of the individual solutions for the vector u> onto the xy, xz, and yz planes 
based on the data from Table 2. Also shown in the figure are the components of our Mean 
1 solution. 

DISCUSSION 

The effects of the actual motions of stars in the differences between the catalogs under 
consideration are ruled out. Therefore, the causes of the rotation c o z found are the following: 

(1) Inaccurate realization of the ICRS or, in other words, residual rotation of the Hippar- 
cos catalog with respect to the ICRF (Ma et al. 1998), and, in this case, the rotation found 
is of a “technical” nature; 

(2) Residual rotation of the ICRF itself with respect to the extragalactic reference frame. 
The ICRF is based on ground-based VLBI observations of extragalactic radio sources. It may 
be assumed that, in this case, the rotation w z found is precessional in nature and depends 
on the accuracy of the constants used. Given the form of the proper-motion differences and 
the signs on the right-hand sides of Eqs. (1)—(2), we can write 

LO y = —A pi sine, (4) 

uj z = Api cose — A E, (5) 

where A pi is the correction to the adopted constant of lunisolar precession in longitude, A E 
is the sum of the corrections to the rate of planetary precession and the motion of the zero 
point of right ascensions, and e is the inclination of the ecliptic to the equator. Using only 
Eq. (5) and setting A E = 0, we obtain from our Mean 1 solution 

Api = —0.38 ± 0.10 mas yr -1 . (6) 

This value agrees with the result of our previous work (Bobylev 2004), where we found 
A pi = —0.42 ± 0.10 mas yr -1 from a kinematic analysis of the proper motions for distant 
Hipparcos stars. In our opinion, a comparison of the residual rotations around the ecliptic- 
axis found by two independent methods is most justifiable. 

Since Ma et al. (1998) adopted the correction Api = —2.84 ± 0.04 mas vr' -1 to the 
IAU (1976) constant of lunisolar precession in longitude, it may be assumed that Api was 
underestimated. Since the ICRS was constructed with precisely this value, solution (6) yields 
an “addition” to this correction. In this case, the correction to the IAU (1976) constant of 
lunisolar precession in longitude is 

A pi = —3.22 ± 0.11 mas yr -1 . (7) 
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Solution (7) is in satisfactory agreement with the most recent results from the laser ranging 
of the Moon, Api = —3.02 ±0.03 mas yr -1 (Chapront et al. 2002), and with the results of an 
analysis of radio interferometric observations, Api = —3.011 ± 0.003 mas yr -1 (Fukushima 
2003). On the other hand, if we assume that the lunisolar precession has no effect on the 
construction of ICRS and ICRF, then A E ^ 0. In this case, we find from Eq. (5) that 

A E = ±0.38 ± 0.10 mas yr -1 . (8) 

The residual rotation of the Hipparcos catalog that we found is small. An analysis of the 
radio-interferometric observations (Boboltz et al. 2003) performed 9.69 years after (epoch 
2000.94) the construction of ICRS (epoch 1991.25) revealed no significant effect in the coor¬ 
dinates of radio stars. However, the rotation components obtained by Boboltz et al. (2003) 
from the coordinate differences for 18 radio stars (Kovalevsky et al. 1997) are of considerable 
interest: 


An = ±0.2 ± 2.9 mas, 

( 9 ) 

e 0y = ±1.9 ± 3.2 mas, 

(10) 

£q z = —2.3 ± 2.8 mas. 

(11) 


We changed the signs of the components of the vector e inferred by Boboltz et al. (2003) 
in order to have Hipparcos catalog differences comparable with those analyzed. To compare 
results (9)—(11), for example, with the Mean 1 solution, the quantities £ox, y ,z must be divided 
by the epoch difference. Thus, the difference e 0z , which is determined with the smallest 
random error, is largest among quantities (9)—(11), thereby confirming our u> z value. In 
general, quantities (9)—(11) agree well with the Mean 1 solution. 

Our analysis of the available individual sources for controlling the ICRS inertiality shows 
that they are of little use for analyzing future projects, such as GAIA, SIM, etc., in which a 
microacrsecond accuracy is expected to be reached (Kovalevsky et al. 1999). The idea that 
the inertiality must be controlled using quasar observations from a spacecraft directly during 
its flight (Metz and Geffert 2004) seems most promising. 

CONCLUSIONS 

We have confirmed that the error in referring the ICRS to the inertial reference frame is 
very small and does not exceed ±0.25 mas yr^ 1 (along the three axes). 

We showed that the equatorial component a j z = —0.35 ± 0.09 mas yr -1 of the vector of 
residual rotation of the ICRS with respect to the inertial reference frame differs significantly 
from zero. This confirms the result of a kinematic analysis of the proper motions for stars of 
the ICRS catalogs (Bobylev 2004). 
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fable I. Components of Ihe vector u >—«; x . ai y , and (in mas yr x ) that we determined from the VLA—Hipparcos 
inferences using data Irom Bobollz et ul. (2003). 


Number of stars 

ob 

1 

^Jjy 

Jj ~ 

15 (all) 

±2.06 

-0.97 ± 0.02 

-1.59 ±0.73 

±0.18 ±0.63 

12 

±1.03 

-0.42 ± 0.56 

-0.51 ± 0.64 

±0.20 ±0.57 


Tlie radio stars HD 50896 N, KQ Pup, and RS CVn were rejected. 


lame z'.\ tqe.UioiVj.ywuvpvMviiKyOVv,'. 
objects— uj x , ujy, uj z (mas yr -1 ). 





P P P 

1 y* 1 z 

AT, 

-^area 

<4r 

Lily 


NPM2 

16.00/4.62/8.16 

3519 

347 

-0.11 ±0.20 

-0.19 ±0.20 

-0.75 ±0.28 

SPM2 

22.15/9.43/28.44 

9356 

156 

±0.10 ±0.17 

±0.48 ±0.14 

—0.17 ± 0.15 

NPMI 

10.24/4.62/16.00 

2616 

899 

-0.76 ±0.25 

±0.17 ±0.20 

-0.85 ±0.20 

SPM1 

44.44/5.71/130.61 

4067 

63 

±0.44 ±0.12 

±0.71 ±0.18 

-0.30 ±0.07 

PUL2 

2.90/1.28/3.63 

1004 

147 

-0.98 ± 0.47 

-0.03 ±0.38 

-1.66 ±0.42 

Kiev(GPMl) 

1.0/1.0/1.0 

415 

154 

-0.27 ±0.80 

±0.15 ±0.60 

-1.07 ±0.80 

Potsdam 

2.37/0.74/2.78 

256 

24 

±0.22 ± 0.52 

±0.43 ± 0.50 

±0.13 ±0.48 

Bonn 

5.54/2.96/5.88 

88 

13 

±0.16 ±0.34 

-0.32 ±0.25 

±0.17 ±0.33 

VLB] 

7.11/2.74/7.11 

12 


-0.16 ±0.30 

-0.17 ±0.26 

-0.33 ± 0.30 

VLA±PT 

2.04/0.45/1.97 

12 


-0.42 ± 0.56 

-0.51 ± 0.64 

±0.20 ±0.57 

HST 

0.08/0.08/0.05 

46 


-1.60 ±2.87 

-1.92 ± 1.54 

±2.26 ± 3.42 

EOP 

8.16/2.36/- 



-0.93 ±0.28 

-0.32 ± 0.28 

- 

Mean 1 




±0.04 ±0.15 

±0.18 ±0.12 

-0.35 ±0.09 

Mean 2 




-0.19 ±0.16 

±0.08 ±0.11 

-0.43 ±0.15 


Note: N* is the number of common comparison stars, and iVarea is the number of comparison areas common to the photographic 
catalogs. Mean 2 was calculated without SPM1. 
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